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I. INTRODUCTION 



Although the case for dark matter (DMHs now quite convincing (for reviews summariz- 
ing the evidence for DM, see, for example, 0]), very little is yet known about its identity 
The mass or masses of dark matter particles remains essentially unconstrained. Aside from 
limits from direct detection experiments [3|-0] and from observations of the bullet cluster 
[§[ , we have essentially no data on DM interactions with itself or with the Standard Model 
(SM). Meanwhile, the list of possible DM candidates is extensive and includes the light- 
est supersymmetric particle, axions, sterile neutrinos, and Kaluza-Klein DM; an extensive 
review of traditional DM candidates is given in j9[. 

Perhaps most importantly, we do not even know how many species of dark matter there 
are, and, even if DM is comprised of just one species of particle, we do not know whether or 
not that particle is just the lightest in some "dark sector" . In recent years, we have seen an 
increasing number of models with more complex DM scenarios, such as two-component DM 
To|. multicomponent DM 11], exciting DM [l2[, inelastic DM [13], and related models 14 



Given this state of affairs, it makes sense to consider models with multiple DM components 
or novel DM-SM interactions. 

At the same time, the SM leaves unanswered many questions regarding the quark and 
lepton flavors. Although in the last few years the precision with which the elements of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix are known has improved significantly (lo| . 
there remains a need for better understanding of the triplication of the fermion generations 
in the SM, or how the three fermion generations could be related through some physics 
beyond the SM (BSM). (For reviews of quark flavor physics, see (l^-18|.) 



Here, we consider the possibility that BSM physics which controls interactions between 
DM and the SM fields and the new physics which explains the flavor structure of the SM 
may be related. We assume that DM belongs to a dark sector which contains at least two 
"flavors" of particles, / and /', both of which we take to be fermionic. 1 We take rrif < rrifi 
and assume that DM is comprised of either / or some mixture of / and /'. We will consider 
a wide range of DM masses, m/ < TeV. 

Given the lack of information on dark matter interactions, and the wide range of flavor 
models available, we find it appropriate to consider SM-DM interactions from a phenomeno- 
logical point of view, instead of trying to incorporate DM into a specific model of flavor. For 
this work, we concentrate on flavor interactions which involve the SM quark fields; as an 
illustration here, we will concentrate mainly on d and s quarks. We consider two example 
interactions: 1) the case where both the SM quarks and the DM fields interact through 
purely vector couplings, and 2) where both interact through purely right-handed interac- 
tions. Generically, one would expect both flavor-conserving and flavor- changing interactions; 
we include both possibilities here. 

The layout of this paper is as follows. First, in Sec. [TTJ we introduce the idea of flavored 
DM and specify our notation and assumptions. In Sec. IHfl we review the current constraints 
on flavor interactions involving DM and d and s quarks from low-energy measurements and 
collider experiments and explore implications of the relic density on the interactions of 
flavored DM with the SM. In Sec. HVt we consider limits from direct-detection experiments, 
taking into account the possibility that the dark sector may contain more than one long-lived 



1 For a model of DM which has implications for flavor physics, but where the DM itself does not carry 
flavor, see [19|. 
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component. Next, we present two toy models of flavor gauge interactions in Sec. |Vj The 
relevance of TeV-scale flavor and DM interactions to LHC is explored in Sec. I VI I Finally, 
in Sec. IVHt we conclude. 



II. NOTATION 



In this section, we will specify the interactions which we will consider among s and d 
quarks and our dark sector particles / and /'. For low-energy observables, we will primarily 
be interested in effective operators of the form 

-J^Olna, = ^(/ml*/»)(?.l>k) (1) 

which give interactions between the SM quarks and the dark sector. We will also occasionally 
consider the four-quark operators 

f9 rig 

-jfOloca- = -^(Qa^q b )(q c T,q d ). (2) 

In these operators, the indices m and n indicate whether we are talking about / or /', while 
a, b, c, and d on the quark fields q indicate whether the quark flavor is s or d. We will assume 
that these flavor interactions are mediated by heavy gauge bosons (which we will denote 
generically as Z'), and, thus, we will confine our discussion to the example cases of purely 
vector (g = V) interactions with r M = 7^ and purely right-handed (g = R) interactions with 

= 7 a»(i + 7 5 )/2. 

Each of these operators is multiplied by a coefficient C 9 mnah j 'A 2 (or C 9 abcd / 'A 2 ) where A 
is taken to be some high new physics scale. As we assume that the SM and DM share 
a common flavor interaction, we expect these operators to have similar scales. However, 
the operator coefficients may also contain small mixing angles; in the SM weak interaction, 
these small mixing angles cause the effective scales between different four-quark operators 
to differ by more than 2 orders of magnitude. Here, we will keep all of our results in terms 
of C 9 mnab /A* and C\J k\ 

For some parts of this analysis, we will have to also consider possible interactions involving 
other quark fields. In the case of vector interactions, we must include both right-handed 
quarks (which are singlets under the SM SU(2) weak interaction), as well as the left-handed 
SU (2) doublets; this necessarily requires that we also consider up-type quarks, weighted by 
the appropriate angles of the CKM matrix. Additionally, the interactions we have specified 
above, if taken in isolation or without careful arrangement of quantum numbers, lead to 
anomalies in triangle diagrams involving the Z's and the SM gauge bosons; however, as 
the choice of quantum numbers for the SM and dark sector fields needed to cancel these 
anomalies is not unique, we will not consider these additional interactions throughout most 
of this paper. 

Finally, among the operators included in Eq. ([2]), there exist some {O g sdsd and their 
Hermitian conjugates) which change strangeness by two units, and, thus, can contribute 
to K° — K° mixing. Constraints on the effective new physics scale for these operators are 
O(10 3 TeV) [181 ] . As we are interested in effects which may be observed at LHC or DM direct 
detection experiments, our analysis will only be applicable to flavor models which do not 
have tree-level contributions to K° — K° mixing. In Sec. |Vj we will present two toy models, 
one with right-handed couplings and one with vector couplings, which are anomaly-free and 
which do not contribute to K° — K° mixing at tree level. 
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III. CONSTRAINTS ON FLAVORED DARK MATTER 



In this section, we will review some of the relevant measurements which constrain BSM 
flavor interactions involving d and s quarks with each other and with the dark sector particles 
/ and /'. We begin with constraints obtainable from low-energy observables. We first 
consider the case where the / (and possibly also the /') is very light, rrif < 180 MeV. In 
this case, we would expect to have the decay K + — >■ n + ff (and possibly decays to final 
states containi ng f or /' as well). Thus, considering the branching fraction of K + to a ir 
plus neutrinos [l5|] 

Br(K + ->■ vr+z/z/) = 1.7 ± 1.1 x KT 10 (3) 

and taking the 2cr bound on this measurement as a limit on the branching fraction to //, 
and taking the ratio of this branching fraction to that of the SM process K + — > e + wK°, we 
obtain, for purely right-handed interactions, 

\r R i i 

A 2 (47 TeV) 2 1 ' 

for rrif (and possibly also rrif) << 180 MeV. We can also consider the case where / is 
very light, but my is somewhat heavier, by comparing to the SM process K + — > fi + u7i°, 
obtaining 

\r R I i 

\^ff'sd\ - I 



< 1 (5) 

A 2 ~ (42 TeV) 2 1 ; 

for rrif ~ m M . (Here, we assume that the limit on K + — > 7i + ff is not substantially degraded 
when one of the final-state dark sector particles acquires a mass of 0(100 MeV).) 

We see from these bounds that the case of very light / is very strongly constrained. 
(Although in both cases here we assumed that / was very light, rrif << 180 MeV, we can 
infer from the strength of these constraints that we would still obtain significant bounds 
on the new physics scale for all cases in which K + decay was not strongly phase-space 
suppressed.) The limits derived here are only applicable to the operators R nsd as shown 
above (and, unless m = n = f, only to the case where both / and /' are light). As 
interaction and mass eigenstates need not coincide, we may be tempted to interpret these 
results as a tentative order-of-magnitude estimate of expected bounds on the new physics 
scale for the operators R ndd and R nss as well; however, small mixing angles (< 0.1) can 
easily invalidate this interpretation. 

Our results are similar for the case of vector interactions, 

\nv I i 

mnsd I ( f\\ 

A 2 ~ O(80 TeV) 2 1 ' 

for the case where both final-state dark sector particles are light, and 

\C V I 1 

mnsd\ fn\ 

A 2 ~ O(70 TeV) 2 1 ' 

for the case where one is light and the other has a mass of O(100 MeV). In this case, 
however, because vector interactions with quarks necessarily involve the left-handed quark 
doublets, and because the interaction eigenstates in the left-handed doublets differ from the 
mass eigenstates by a rotation via the Cabibbo angle 9c, the operators 0^ nab necessarily 
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involve some significant quark mixing. Since we have no good reason to believe that the 

interaction eigenstates would be aligned with the down-type quarks, it therefore may be 

\c v I 

more compelling in this case to interpret the bounds on — v j^ L - as a tentative bound on the 
a = b = d, s cases as well. 

We note that one can also obtain limits on the right-handed operators for the case a = 
d from limits on supernova cooling, assuming that rrif < few x 10 MeV. We can apply 



the limit on the branching fraction Br(7r° — > vv) < 3.2 x 10 -13 [15|, |20j to the case of a 
light /. For rrif m few x 10 MeV, we obtain limits on the new physics scale of order ~ 1 
TeV. However, as the decay 7r° — > ff is helicity-suppressed, this limit quickly degrades with 
decreasing rrif. Similar statements apply in the case that one or both final-state particles is 
an /'. 

We now move on the the case of heavier /. We will first consider constraints which can 
be obtained from K° — K° mixing. As we limit ourselves to models which have no tree-level 
contribution to K° — K° mixing, we will only concern ourselves with contributions from the 
diagram shown in Fig. [TJ (In the case of vector interactions, there are also loop diagrams 
containing up-type quarks which can contribute to K° — K° mixing. We will briefly discuss 
these contributions when we consider toy models in Sec. |Vj) For simplicity, we will only 
consider the case where rrif,mfi » m K /2. 




FIG. 1: K° — K° mixing contribution from internal dark fermion loop. 

First, we will consider operators which change flavor on both the SM and dark sector 
fields. For a right-handed interaction, we have 

-jjr°rf*. = -^U'nl'mdRl.sn) (8) 
while for a vector interaction, we have 

= % £ (/V/)(rf7^). (9) 

First, we note that these interactions, taken in isolation, do not contribute to the diagram 
in Fig. (TJ the / created in one of the vertices in the loop would have to transform into an /' 
before it is destroyed at the other vertex. However, we note that, in general, we would not 
expect the interaction eigenstates in the dark sector to necessarily coincide with the mass 
eigenstates / and /'. We can introduce interaction eigenstates fi and /2, related to / and 

f r)-UZ)ti) 



where \a\ = 1 and, in the case of the right-handed interaction, the fields in Eq. ffTOj) are 

taken to be the right-handed components only. (Here, we assume that / 1)2 do not mix with 
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any other states.) We can now place limits on combinations of the coefficients Cj 2 j ids / A 2 , 
the parameters a and ft which control the mixing, and the mass splitting, 5 = my — m/. 
In the case of right-handed interactions, the diagram in Fig. [1] is finite; we obtain, for the 
mass difference between K L and K$, 



A, 



m K 



Af z K m K a* 2 ft 2 



r<R 

L/ f 2 fid s 
A 2 



(4tt) 5 



(11) 



where fx ~ 160 MeV is the kaon decay constant and A is a constant which depends on the 
relative size of mj and my] for m/ w m/i, A « 8/9, while, for m/ << my, A 2/3. We can 
then compare this to the experimental value A mK = 3.48 x 10~ 15 GeV 15); for A = 1 TeV, 
we obtain upper bounds on \Cf 2 j ids a* ft8\ of approximately 7 — 8 GeV. 

We get a similar order-of-magnitude expression for the case of vector interactions, 



A, 



m K 



m 3 K 



3 JK (m s + m d ) 2 



a 



*2 n2 



2 



A 2 



(47T) 



r5 2 ln 



A 2 



m. 



(12) 



where the logarithmic behavior comes from the running of Cj sds between the new physics 
scale and my induced by the diagram in Fig. [TJ For A = 1 TeV, this gives an upper bound 
on \CYf lds a*ft8\ of ~ 1 GeV. Here, we have taken the logarithmic factor in the vector case 
to be of order unity, but we note that the upper bound on \Cj 2 f lds a* ftS\ can be strengthened 
considerably if / and /' are very light. Thus, we see that fairly small mass splittings are 
phenomenologically interesting in these scenarios. 

We also briefly consider the operators which do not change dark sector flavor. For both the 
right-handed and vector interactions, the diagram in Fig. [T]will be quadratically divergent. 
If we regulate this divergence using dimensional regularization, we find that the contribution 
to A mK from the operator for the right-handed case, 



U ffds ^R 



C R 



(13) 



will contain factors of mj instead of 5 and, thus, this operator will be more strongly con- 
strained than the operators which change dark sector flavor as long as my is substantially 
larger than 5. The analogous vector operator, on the other hand, will contain factors of tuk 
instead of S or m/; the scale for this operator will be constrained to be at least O(TeV). 

We also mention that the limits derived from A mK above depend only on the magnitude of 
^%hds a * ft, n °t its phase. However, the CP-violation parameter is sensitive to a complex 
phase in the kaon mixing matrix, which depends on Cj 2 j ids a* ft . Depending on the choice 
of these phases, the upper bound on \Cj 2 j ids a* ft5\ could be strengthened by more than an 
order of magnitude for a given value of A. 



We briefly mention relevant collider bounds. CDF [2l| has directly looked for heavy 
neutral gauge bosons, Z', decaying to jets. They have excluded the mass range of 320 — 740 
GeV, assuming SM couplings. (For limits on lighter Z' gauge bosons, see results from UA2 



23 



However, these limits can be easily evaded by making the couplings between the 
quarks and the Z' slightly smaller than those in the SM. Additionally, indirect limits on 
specific flavor models from fits to collider data were given in [IT 
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Fin ally , we consider constraints obtained from the observed DM density of tldark — 0.228± 
0.013 [25(1 . If DM consists of only one species, obtaining the correct relic density requires 
a velocity-averaged annihilation cross-section at freezeout < av r > of ~ 3 x 10 _26 cm 3 /s, 
where v r is the relative speed of DM particles, with a mild dependence onmj. As long as 
rrif > O(GeV), each of our purely right-handed operators will contribute a term 

\fiR 12 3 m 2 

to the annihilation cross-section f m f n — > q a (jb- (Here, we have neglected velocity-dependent 
terms.) This gives a value for the sum of the squares of the operator coefficients, 



E 



\n R I 2 1 /TpV\ 2 

niiudA 1 / 1 4 1 (15) 



A 4 (2.6 TeV) 4 \m f J 



We see that the new physics scale decreases with decreasing mf. We note, however, that, 
if we allow / to comprise only a fraction of DM (and thus have a smaller relic density), we 
can allow larger values for < av r >, and thus lower values for the new physics scale. Thus, 
if we have only one operator with a nonzero coefficient, say, Off M , then 2.5 TeV should be 
regarded as an approximate upper bound on the new physics scale. (Note that all of our 
other constraints place lower bounds on the scale of new physics.) Of course, this upper 
bound on the scale of new physics can be raised by having additional terms with nonzero 
coefficients in the sum shown in Eq. ( 1151) . so this should not be taken as a rigorous upper 
bound on the NP scale. 
We get a similar term 



v \C V h \ 2 6rn 2 f 
< t7v r > v mnab = ^f^^ (16) 

for the operators with purely vector interactions. Here, we have included a factor of 2 in 
the cross-section, as gauge invariance demands that we include both the upper and lower 
components of the left-handed quark doublets. This gives 



E 



\n v I 2 1 /TpV\ 2 

nnuilA 1 / 1 v 1 (17) 



A 4 (5.2 TeV) 4 \m f J 



We note that our relic density calculations assume the simplest scenario for DM self- 
annihilation // — > qq. If / and /' have masses which differ by only a few percent or less, the 



more complicated coannihilation calculation [26J is relevant; coannihilations can significantly 



reduce the relic density of /, and thus loosen the limits in (ITS]) and (fTTj) . if the annihilation 
cross-section for /'/' is substantially larger than that for //. Additionally, this calculation 
also neglects the possibility of other annihilation channels for //, which can substantially 
raise the new physics scale. 



IV. SIGNATURES IN DARK MATTER DIRECT DETECTION EXPERIMENTS 

We now discuss the prospects for direct detection of flavored DM. We will distinguish 
between two experimentally distinct cases, depending on the mass splitting 5. In the first 
case, 5 > few x 100 keV — 1 MeV. (The range of masses given here reflects the difference 
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between the right-handed and vector operators.) In this case, / is the only long-lived particle 
in the dark sector, and it interacts elastically with nucleons, fN — > fN. In the second case, 
as we will show below, it is possible that the heavier state /' can also be long-lived and 
form a sizable component of DM, thus allowing exothermic down-scattering of the form 
f'N —> fN. (Additionally, for 5 < few x 100 keV, up-scattering of the form fN —> f'N 
becomes possible, thus making flavored DM a possible example of inelastic Dark Matter 
introduced to possibly explain the apparent conflict between DAMA [27| and other 
experiments. Although this scenario is significantly constrained 0, l28l. [29| , it remains a 
viable possibility for explaining the DAMA results if 5 ~ 0(200 keV) [30J. 2 ) 

First, we will discuss the case where 5 is sufficiently large that /' is not long-lived, and 
DM can interact in direct detection experiments only through the elastic reaction fN — > 
fN. In this case, we can directly apply the constraints from existing DM experiments; 
here, we will only consider spin-independent contributions to the cross-section, as they 
are significantly more strongly constrained than spin-dependent contributions. As direct 
detection experiments search for interactions between DM and nucleons, we are interested 
in operators which contain d quarks (and u quarks, in the case of vector interactions). For the 
purely right-handed interactions, the measured spin-independent DM-nucleon cross-section 
takes the form 

asi ~ ~l^T^ M - d a* ' (18) 

where M re d is the reduced mass of the DM-nucleon system and A and Z are the atomic mass 
number and atomic number of the target nucleus, respectively. The dependence on A and 
Z takes into account the fact that limits from DM direct-detection experiments assume that 
the cross-sections for DM scattering on protons and neutrons are equal; / interacts only 
with d quarks, and not u quarks; thus, the cross-section for a neutron is four times that of 
a proton. For the purely vector interactions we have 

*8i = -^-M^ (19) 

where the factor of 9 includes the contributions from both u and d quarks. (As / interacts 
identically with u and d quarks, the cross-sections here are the same for protons and neutrons; 
we ignore corrections due to the Cabibbo angle.) 

We can compare these expressions to the cross-section limits from direct detection experi- 
ments; here we will assume that all of DM is comprised of /. We consider three specific values 
for rrif. For the CoGeNT signal region, rrif « 10 GeV [31], we take a si ~ 5 x 10 -41 cm 2 
(although this is in conflict with [5|, |6j). For mj ^ 1 TeV and the range of DM masses 
where these limits are strongest, rrif m 0(lO's of GeV), we use the results from XENON100 
who report a spin- independent cross-section upper bound for these mass ranges of 
&SI ^ 8 x 10~ 44 cm 2 and asi < 7 x 10 -45 cm 2 , respectively. Our results are shown in 
TableE 



2 



We note, however, that arranging for up-scattering to occur while suppressing the elastic interaction 



fN — > fN presumably requires significant fine-tuning of mixing angles; [30 ] utilize inelastic scattering 
cross-sections which are 10 orders of magnitude larger than the elastic scattering cross-sections ruled out 
by 0] to address the DAMA results. This does not preclude, however, the existence of inelastic scattering 
with much smaller cross-sections, which could require less fine-tuning of mixing angles. 



s 



m f 


(|Cf f(M |/A 2 )-VVTeV 


(I^IM 2 )- 1/2 /TeV 


~ 10 GeV 


~ 0.7 


~ 2 


few x 10 GeV 


> 7 


>19 


~ 1 TeV 


> 4 


>H 



TABLE I: Results for the new physics scale from DM direct detection experiments. The first line 
refers to NP scale corresponding to the CoGeNT signal region [3fl, while the last two are lower 
bounds on the NP scale obtained from the results of XENON100 lal. 



It should be noted that these direct detection limits on the lower bound for the new physics 
scale disagree with the upper bounds from the relic density calculation unless / interacts 
with more than just d (and, in the case of the vector operators, one family of up-type) quarks. 
Given the significant tension between the direct detection limits and the relic density bounds, 
it may be necessary to solve this problem by, for example, introducing annihilation channels 
of DM particles into leptons [31], introducing additional interactions and/or particles to 
the dark sector, or positing mechanisms which would give / a non-thermal cross-section. As 
these issues are not specific to flavored DM, we will not try to address them here; we will 
instead take the attitude that explaining the relic density will require the construction of a 
specific model. (We also note that if we relax the assumption that / comprises all DM, the 
direct detection constraints are also loosened.) 

We now address the question of when the excited state /' can be suitably long-lived 
to form a significant component of DM. Here, we make the conservative assumption that 
the couplings of //' to dd and ss are not significantly smaller than those to ds and sd. 
(If couplings which change SM flavor dominate, /' decays will be suppressed by factors 
greater than those shown below for the flavor-diagonal case; this may happen, for example, 
in the toy models presented in Sec. |V] if the DM interaction eigenstates are closely aligned 
with the mass eigenstates.) For the case of right-handed interactions, /' can decay at tree- 
level as /' — > fn° if 5 > m n ; tree-level decay for the vector coupling case can occur via 
/'—>•/+ jets for somewhat higher 5. In both of these cases, /' will not be long-lived unless 
the scale of new physics is extremely high; thus, we consider the cases of substantially smaller 
8 and decays involving neutrinos and photons. For the case of right-handed interactions, we 
consider /' — > fvv and decays containing one or two photons. We will find that /' decays 
are substantially more weakly constrained in the case of vector couplings; in this case, we 
will consider the decays /' — » ^777 and /' — > fe + e~ as well. We note that, while we are 
interested in the case where /' lives long enough to comprise a significant fraction of DM 
today, there also exist constraints on DM particles with lifetimes substantially less than the 
age of the universe but greater than 0(1 second); see 33[ and references therein. 

For the case of right-handed interactions, the strongest constraint we obtain is by con- 
sidering the /' /7 diagram shown in Fig. [2j This diagram is superficially logarithmically 
divergent; however, if we consider this diagram from an effective-operator point of view, 
gauge invariance requires that the decay occur via an effective transition magnetic moment 
operator, fcr^f'F^. As this operator contains an explicit factor of the photon momentum, 
the number of factors of loop momenta that we integrate over must be reduced by one; this 
renders the diagram finite. We will assume that the integrals over the loop momenta can be 
reasonably estimated by cutting them off at our new physics scale A; for simplicity, we take 
the C% nab to be 0(1). 

The diagram in Fig. [2] contains a sub-diagram which couples two Z' bosons to a photon 
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FIG. 2: Diagram contributing to /' — > f*y in the case of right-handed couplings to the Z' . 



via a quark loop. We take this sub-diagram to be similar to the fermion loop contribution 
to the SM effective ZZj vertex jlij], which depends very weakly on the fermion mass in the 
scenario of large gauge boson momentum, as is relevant here. However, this sub-diagram is 
anomalous (as will be discussed in more detail in Sec. [V}. Therefore, when we sum over all 
possible fermions in the loop, terms which are independent of fermion mass will cancel, and 
we must look at mass- dependent terms, which, roughly speaking, will contribute a factor of 



^M 2 

arrive at our estimate for the decay of /' to fj: 



rriq/A to the decay amplitude, where m q is the mass of the fermion in the loop. We thus 



(4tt) 8 \A/ A 4 



We then utilize the results of 35], who use the Milky Way 7 ray line search to constrain 
dark matter decays containing monoenergetic photons. For the values of 5 relevant here, 
they exclude such decays if V > 10 _49 mj. Here, we take A = 1 TeV and m q = m s 3 ; this 
latter choice is valid if the flavor charges of the SM are arranged to cancel the anomalies; if, 
instead, new (heavy) particles are added, substantially stronger constraints will be obtained. 
For these values, we obtain 

6 < (1200 keV) (—) (21) 

as an order-of-magnitude upper bound on the values of 5 which will allow /' to be long- 
lived. Although this is an approximate limit, we would like to point out that considering /' 
decaying via a virtual pion will give values for 5 of the same order of magnitude. 

For vector interactions, the /' lifetime depends strongly on whether or not the decay 
channel /' — » fe + e~ is kinematically allowed. For 5 < 2m e , only the decay channels /' — > 
fuu and /' — > f + photons are allowed. We considered many diagrams that could contribute 
to these decays. The Z' does not mix into an on-shell photon, so /' — > does not occur at 
one-loop order. A quark loop connected to three vector bosons vanishes identically, which 
eliminates /' — > f^py at one-loop order as well as /' — > /7 via the two-loop diagram shown 



We take the value of ms(l TeV) = 47±j4 from [36|; we note that the error on this number can probably 



be reduced using the recent results of [37[, but we emphasize that here we are interested in an order-of- 
magnitude result for 6. 
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in Fig. [21 /' — > /77 via a virtual ir° does not occur for a purely vector interaction, and 
/' — > /77 via a virtual p is forbidden by charge-conjugation invariance. The largest nonzero 
contributions to /' decay via Z's are /' — > fvi> via Z' — Z mixing and contributions to 
/' — > /777 via a quark loop or virtual p and/or tt° mesons; example diagrams contributing 
to these processes are shown in Fig. [3J Both of these processes will be strongly suppressed; 
/' — > jvv will be suppressed by factors of both the Z' and Z masses, while the rate for 
/' — > /777 is suppressed by many factors of the small photon momenta. As in the case 
of right-handed couplings, if we do not add additional fermion fields, anomaly cancellation 
requires that coefficients of the operators with first and second generation quarks have equal 
magnitude and opposite sign, which renders the diagram in Fig. [3]^a) finite. If we take 5=1 
MeV and insist that /' have a lifetime comparable to the age of the universe, this diagram 
gives limits on the new physics scale weaker than O(GeV). (We note that it may be possible 
to slightly improve these limits using observations of dark matter halos 38].) We also obtain 
order-of magnitude constraints on the new physics scale using diagrams for /' — > ^777 such 
as those in Fig. [3] (b) and using the limits on dark matter decays involving photons from 
[35I ]; for 5=1 MeV and rrif ~ 100 MeV, the limit on the new physics scale is no better than 
O(10 GeV); this limit weakens with growing m/. 






FIG. 3: Diagrams contributing to /' decay in the case of vector couplings to the Z'. a) /' 
via Z' — Z mixing, b) Example diagrams contributing to /' — > /777- 



Although the limits from these diagrams are very weak, we note that it is possible to 
have the /' decay in this scenario for 5 < 2m e and for interesting values of the new physics 
scale if we make some assumptions about the scalar sector of the model. We will mention 
this briefly in Sec. IVl when we discuss toy models. 

Finally, we mention for the case of vector interactions that if 5 > 2m e , f can decay much 
faster via a diagram similar to that in Fig. [3] (a) but with the SM Z replaced by a photon 
and with the neutrinos replaced by an e + e~ pair. 

Thus, we conclude that for 5 < few x 100 keV (for right-handed interactions) or 5 < 
1 MeV (for vector interactions), /' can possibly be long-lived. This opens up the possibility 
that /' could be discovered via its distinctive down-scattering signatures at direct-detection 
experiments; for works where such signatures have been considered, see (H-lii. Although 



these signatures would be useful only for very small 5, if observed, they would provide strong 
evidence for multicomponent dark matter. 
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V. TOY MODELS 



Here, we present two toy models which contain, respectively, the purely right-handed and 
purely vector couplings described above. We show that they are anomaly-free and do not 
contribute significantly to K° — K° mixing. In both cases, we assume that quarks from the 
first and second generations transform as a doublet under a gauged SU(2) p flavor symmetry, 
under which all other SM fields are singlets. In the first case, only the right-handed down- 
type quarks transform under the SU(2) p flavor symmetry, while in the second case all quarks 
from the first and second generations transform under the symmetry. In both scenarios, the 
/ and /' also transform as a doublet under the same SU(2)p symmetry. We now describe 
the particle content and couplings in each of these toy models. 



A. Right-handed couplings 

In this model, the right-handed down-type quarks form a doublet under SU{2)p: 

dm 



IR2 



Dn =\J L ) (22) 



while all other SM fields are SU(2)p singlets. Here, we perform a rotation such that d\ 
and di are aligned with the right-handed components of the mass eigenstates dn and sr, 
respectively; we neglect possible mixing with 6r. Similarly, we have 



'R 



fm 

fR2 



(23) 



where we do not assume that fm and /r2 are necessarily closely aligned with the mass 
eigenstates /# and f' R . The three gauge bosons of SU(2)p are labeled Z[ for i = 1,2,3; 
we assume that these gauge bosons obtain their mass through a scalar SU(2)p doublet if 
acquiring a vacuum expectation value which, for some SU(2) rotation, can be written in the 
form 

* ( 30 " (24) 

We must consider the possibility of anomalies arising from triangle diagrams in this 
model. For a given triangle diagram with external gauge bosons A^, A b , A° p , the anomaly is 
proportional to 

Tr[(-l) n T a {T b ,T c }) (25) 

where the T % are the generators corresponding to each of the gauge bosons, n = 0(1) for 
left-handed (right-handed) fermions, and the trace is over all fermions that can run in the 
loop. A triangle diagram with exactly one Z' boson does not give a triangle anomaly, as 
the trace over a single SU(2) p Pauli sigma matrix in ( 12 5 p gives zero. A diagram with three 
external Z' bosons similarly vanishes. 

Thus, we need only consider those diagrams with two or zero SU{2) p gauge bosons. For 
those diagrams with two SU(2)p gauge bosons, Z[ and Zj, those which contain a graviton 
or a SM SU(3) or SU{2) gauge boson are zero. However, there is a constraint from the 
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diagram with a hypercharge gauge boson; the anomaly for a diagram with an external Z[, 
Z'j and hypercharge gauge boson B is proportional to 

TT[(-l) n Y{r\r^}] (26) 

where Y is hypercharge, and the r a are the SU(2)p Pauli sigma matrices. Setting this to 
zero implies 

5 ij (-l) n ^ = (27) 

SU(2) F 

where the sum is over SU(2)f doublets. With just /, /', and the SM particle content, 
this relation is not satisfied, as the only nonzero term in the sum is Dr, with hypercharge 
Y = —1/3. We can solve this problem by adding either another right-handed SU(2)p 
doublet with Y = 1 or a left-handed SU{2) F doublet with Y = — 1 which is a singlet under 
the other SM interactions. (Here a factor of three arises since quarks carry color.) However, 
we must also make sure we do not produce anomalies via diagrams which contain no Z' 
gauge bosons, as this would spoil the anomaly cancellation of the SM. We can achieve this 
by adding in two additional SU(2)p right-handed singlets, with hypercharge Y = —1. We 
assume that all non-SM fields other than / and /' are sufficiently heavy to have escaped 
current experimental constraints. 

We now consider the possible contributions to K° — K° in this toy model. At tree level, 
the operator Of dsd = s^^ds^^d is not generated, as it does not obey the SU{2) symmetry. 
Additionally, loop diagrams which do not contain SM W bosons or the DM particles /, /' 
will not generate Of dsd . Although one might expect this operator to possibly arise at one 
loop due to SM W ± exchange, no contributions arise beyond those of the SM; as Of dsd 
contains only down-type quarks, any diagram which contains only one W boson must be 
such that the W starts and ends on the same quark line. However, as the Z' does not 
couple to up-type quarks, no one-loop diagram with only one W can be constructed which 
contributes to Of dsd . With two W ± bosons, one recovers the usual SM K° — K° mixing 
contribution. However, when we include the diagram shown in Fig. [TJ we obtain the results 
described in Sec. IIHI (Similar conclusions also apply to the SU{2) F doublet added in to 
cancel anomalies, as discussed above.) 



B. Vector couplings 

In this case, both the right-handed and left-handed quarks transform under the SU(2)p 
symmetry. As the left-handed quarks come in 577(2) weak doublets, we must include both 
the up- type and down- type quarks here. Thus, we define, along with Dr, 

Ur =( UR1 ) (28) 



and 



UR2 



Qli 



Q L =[q l2 ) ( 29 ) 



where the Qu are the left-handed weak quark doublets. Although we again take the inter- 
action eigenstates to be quasi-aligned with the mass eigenstates, there is necessarily some 
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mixing with the third generation, as the weak quark doublets are themselves not precisely 
aligned with the mass eigenstates. 

We now consider anomalies for the case of vector interactions. We again must consider 
the triangle diagram with two SU(2) p gauge bosons Z[ and Z'^ and one hypercharge gauge 
boson. In this case, the constraint is 

8 ij {-l) n mY (30) 

SU(2) F 

where notation is as before, except that m = 1 for Ur and Dr, but m = 2 for Ql, as Ql 
contains both up- and down-type quarks. Y{Ur) = 2/3, Y(Dr) = —1/3, and Y(Ql) = 1/6. 
Thus, the sum in Eq. fl30l) is zero without the addition of any particles beyond the SM. 
Therefore, we need not consider diagrams containing only SM gauge bosons, as the anomalies 
cancel just as they do in the SM. 

We now consider possible contributions to K° — K° mixing. Like in the right-handed 
case, 0^ dsd is not generated at tree level. However, it can be generated at one-loop level 
via the diagram shown in Fig. [T] or by a similar diagram with the dark sector particles 
replaced by quarks, or by the diagram in Fig. HI The diagrams containing quark loops are 
suppressed by two factors of SM quark mixing angles and by quark masses (with the largest 
contributions coming from c quarks) and their contributions to A mx give limits on the new 
physics scale of 0(1 TeV). 




FIG. 4: K — K mixing contribution in the case of vector interactions. 

We found in Section DID that K° - K° mixing implied a small mass difference 5 or a small 
mixing angle between the mass and interaction eigenstates. For vector interactions, the 
SU(2)p symmetry allows a bare mass term mFF without a coupling to <p. Mass splitting 
and mixings can be accomplished, however, by coupling FF to (p through higher- dimensional 
operators, such as Fipip^F. As we would expect these operators to be suppressed by some 
mass scale, it may be reasonable to generate small mass splittings (even for large rrif) and/or 
small mixings in the dark sector. 

If we assume that such higher-dimensional operators give effective couplings of ip to both 
our DM sector and the quark sector, we can have diagrams that allow /' to decay, in addition 
to those studied in Sec. IIVI Here, we consider a diagram similar to that in Fig. [2] but with 
one of the Z' bosons replaced with a scalar; we take the quark loop to contain a c quark. If 
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we take the effective coupling yp of ip to FF to be roughly ~ S/v', and also assume that a 
similar relation holds for the quark mass splittings, such that y c ~ m c /v', we estimate the 
decay width of the /' to be 

~ m f \4:7Tm f ) V (4vr) 8 A 4 J [ ' 

where the term in the first set of parentheses is a phase space factor and the term in 
the second set of parentheses comes from an order-of-magnitude estimate of the amplitude 
squared for this two-loop diagram. For 5 = 1 MeV and v' ~ A ~ TeV, this gives T ~ 10 -50 
GeV; for this mass splitting, the results of jHf exclude an /' — > f'y for which T > 10~ 49 m/; 
thus, for rather small rrif, this diagram starts to probe an interesting region of the new 
physics scale. We find no diagrams which give stronger constraints. 



VI. SIGNATURES AT LHC 

We now consider the prospects for discovering flavored DM at LHC. For the new physics 
scales which we have been considering, 0(few TeV), the effective operator formalism which 
we have been using up until now is no longer applicable. Here, we assume that flavor 
interactions are mediated by a heavy gauge boson, which we will label as Z'. We note 
that the new physics scales which we have been considering include couplings, and that, for 
example, a new physics scale of 3 TeV, which is allowed by K — K mixing and can still be 
accommodated by direct detection experiments, could easily correspond to a 1-TeV Z' . 

For concreteness, we will consider the right-handed toy model described in Section |V] 
This model contains three Z' gauge bosons, all of which we will generically denote Z' , and 
all of which we will assume have a mass of Mz' = 1 TeV. We will take the / and /' to 
have masses much less than Mz>, and, for simplicity, we take the mixing angle between 
interaction and mass eigenstates to be zero. All other fermions which are added into the 
model to cancel anomalies are assumed sufficiently massive that they are not accessible at 
LHC. The SU(2) coupling is taken to be the same as the SM £77(2) coupling g. These 
assumptions correspond to values of the effective operator coefficients 

1 



\r<R i 
\ u ffdd\ 


\r*R i 


\r*R i 

\^f'f'dd\ 


\r*R i 


9 2 „ 


A 2 


A 2 


A 2 


A 2 


' 4 TeV 2 ~ 






\r<R i 

l^ff'sdl 


\r<R i 

l°/7<isl 


g 2 






A 2 


A 2 


2 TeV 2 ~ 



(3 TeV) 2 
1 

(2 TeV) 2 



(32) 



with all other coefficients C^ nab between the dark and SM sectors 0. 

As an example of a possible search channel at LHC, we will consider the case where the Z 1 
is produced in conjunction with a jet, pp — > Z'j 4 , and decays invisibly to produce a monojet 
signature. (For a discussion of m onoj et signatures at LHC, see [44|; for previous work on 
monojets with regards to DM, see [45j.) An invisible decay of the Z' can consist of //, but, 
in the case where /' is adequately long-lived to leave the detector (which we assume here), 



4 Here, we consider all hard subprocesses of the forms q a g — > q^Z' , q a g — > q^Z' , and q a qi, — > Z'g where q a 
and qb can be either s or d. 
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we must also include the final states /'/, //', and /'/'. We use MadGraph/MadEvent [46[ 
to calculate the cross-sections for these processes and for the SM backgrounds pp — > Zj 
(with Z — > vv) and pp — > W ± j (where the W decays leptonically and the charged lepton 
is lost down the beampipe, r\ > 2.5). We greatly reduce these backgrounds by placing a very 
tight cut on the transverse momentum p? of the jet. 

The pt distributions for our Z' signal and the SM background are shown in Fig. for a 
center-of-mass energy ^fs = 14 TeV. Although the two distributions are similar, the signal 
distribution falls off more slowly for large px- Given the eventual expected data set for 
LHC (an integrated luminosity of ~ 100 fb _1 ), the discovery potential for such a Z' would 
be expected to be limited by systematic errors, not by statistics. Although a full study of 
systematics is beyond the scope of this paper, we note that a value of S/B of 10% can be 
obtained by requiring the monojet to be > 440 GeV. This cut has been applied in Fig. 
El giving signal and background cross-sections of 0.047 pb and 0.47 pb, respectively, and, 
for an integrated luminosity of 100 fb~\ S/^B « 22. S/B = 20% can be achieved with a 
Pt cut of 625 GeV. Although this is far from a complete analysis of this signature at LHC, 
these numbers indicate that this search channel merits further study. 




Pj/GeV 

FIG. 5: Monojet px distributions for a 1-TeV Z' signal and the SM background at yfs = 14 TeV. 
A cut has been placed requiring > 440 GeV, which gives S/B = 10%. 

Wc briefly mention a few other possible signatures of this model at LHC. First, we note 
that have studied the potential of discovering an invisibly-decaying Z' produced in 
conjunction with an SM Z which decays to leptons. They conclude that a 5a discovery of a 
1-TeV Z' from a BSM U (1) with gauge coupling of unity could be accomplished with 30 ftV 1 
of data. (See also [48] for a similar study of an invisible Z' produced in conjunction with 
a 7.) Additionally, visible decays of the Z' could be considered. In addition to Z' — >■ jj, 
we have, for example, Z' — > /'/'; for my > few GeV, and m/ << my, an /' with a few 
hundred GeV of energy will decay in the detector, which can give, among other signals, 
Z' — > j'f — > ffjjjj; depending on m Z ', rrif, and rrif, displaced vertices are also possible. 
We do not consider these signals here, but mention that they could be studied in a more 
complete treatment. 
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VII. CONCLUSIONS 



We see that flavored DM is a rich subject. Here, we have examined DM which interacts 
with quarks of the first two generations; we have then placed constraints on these interactions 
using low-energy measurements and direct detection, and considered the implications of the 
relic density for possible flavored DM models. We also see that flavored DM has possible 
signatures at LHC and that it can give inelastic scattering (both up- and down-scattering) 
in direct detection experiments. Throughout this analysis, we have strived to be as model- 
independent as possible. 

We collect some general results in Table (jlT])- Here, we show the approximate NP scale 
probed by each of these observables under the assumption of flavored DM. We must empha- 
size that not all of these results will apply to all models, and, without a specific model, these 
numbers should not be compared to each other. (For example, the results from K + decays 
are only applicable to very light DM.) Additionally, these results can also be significantly al- 
tered in models with small mixing angles. For these reasons, we give only order-of-magnitude 
estimates of the reach for each of these observables. We do not include a number for K — K 
mixing as it depends very strongly on the mass splitting 5. Additionally, we have not done 
a complete study of the signatures of flavor Z's at LHC, but we take a few TeV to be a 
reasonable estimate of the new physics reach for these scenarios. 



Observable 


Approximate NP scale reach 


K + decays 


O(40 - 80 TeV) 


K — K mixing 


(^-dependent 


Relic density 


few TeV 


Direct detection (elastic) 


0(1 - 10 TeV) 


LHC 


few TeV 



TABLE II: Order-of-magnitude estimates of the NP reach for various observables in flavored DM 
scenarios. Note that, without a specific model, these numbers cannot be meaningfully compared 
with each other. For more detailed information, see text. 

Under the assumption that our particle / comprises all of DM and couples to first- 
generation quarks, the constraints from direct detection are quite strong; in particular, in 
the case of vector couplings, only a very light / would possibly be observable at LHC. 
However, considerations of the relic density indicate that additional interactions may be 
necessary; these additional interactions may introduce new signatures at LHC. Thus, it may 
be fruitful to attempt to incorporate flavored DM into a more complete model. If such 
a model contained multiple types of DM, the constraints on couplings to first-generation 
quarks from direct detection could also be loosened, opening up additional parameter space 
accessible at LHC. 

Finally, we note that we have confined ourselves to interactions involving d and s quarks 
(plus u and c quarks where necessary), and have limited ourselves to purely right-handed 
and purely vector couplings. However, we would like to point out that the range of flavor 
interactions which could potentially be applied to DM is immense. One could consider 
scalar interactions, interactions with leptons, and, perhaps most interestingly, interactions 
involving the third family of quarks. This last option in particular could lead to interesting 
signatures in top physics at LHC. We leave these ideas for future work. 
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